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Abstract 

A systematic analysis of the multifragmentation (MF) in fully reconstructed 
events from lA GeV Au, La and Kr collisions with C has been performed. 
This data is used to provide a definitive test of the variable volume version of 
the statistical multifragmentation model (SMM). A single set of SMM param- 



eters directly determined by the data and the semi-empirical mass formula are 
used after the adjustable inverse level density parameter, is determined by 
the fragment distributions. The results from SMM for second stage multiplic- 
ity, size of the biggest fragment and the intermediate mass fragments are in 
excellent agreement with the data. Multifragmentation thresholds have been 
obtained for all three systems using SMM prior to secondary decay. The data 
indicate that both thermal excitation energy E^f^ and the isotope ratio tem- 
perature Thc-dt decrease with increase in system size at the critical point. 
The breakup temperature obtained from SMM also shows the same trend as 
seen in the data. The SMM model is used to study the nature of the MF phase 
transition. The caloric curve for Kr exhibits back-bending ( finite latent heat) 
while the caloric curves for Au and La are consistent with a continuous phase 
transition ( nearly zero latent heat) and the values of the critical exponents r, 
P and 7, both from data and SMM, are close to those for a 'liquid-gas' system 
for Au and La. We conclude that the larger Coulomb expansion energy in Au 
and La reduces the latent heat required for MF and changes the nature of the 
phase transition. Thus the Coulomb energy plays a major role in nuclear MF. 
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I. INTRODUCTION 



During the past decade a large effort has been made to understand the multifragmentation 
(MF) process in heavy ion reactions. Recently a number of review articles have appeared 
describing the details of this phenomenon The Purdue group was the first to suggest 

that nuclear MF might be a critical phenomenon - a second order phase transition occurring 
near a critical point P,^. In an inclusive experiment performed at Fermilab, the spectra of 
mass identified fragments resulting from proton interactions with both Kr and Xe gas targets 
were measured Pfl- The discovery that the yields of fragments with mass Af produced in 
p+Kr and p+Xe reactions, with proton beam energy from 80 to 350 GeV, obeyed a power 
law Y{Af) ~ ^/^) with r ~ 2.5 generated theoretical interest in MF in terms of a 
continuous phase transition. A similar power law was also predicted by Fisher for a 
mass distribution of droplets at a liquid-gas phase transition critical point. Thus the results 
from the Purdue work gave a hint that MF could provide important information about the 
equation of state of nuclear matter |]TT]-|T^. The analysis of the fragment kinetic energy 
spectra suggested that the fragments are emitted from a less than normal density system in 
the decay of a common remnant which is lighter than the target [p|,|r^-|T7|| . An analysis of 
the fragment yields based on a thermal liquid drop model gave a freeze-out temperature of 
~ 5 MeV Similar results were observed when the above study was extended to lower 



energies [|T8|,|T9[. This view was supported by exclusive emulsion MF data ||2^, which were 
analyzed by Campi to show that the conditional moments of the individual fragment events 
exhibited characteristics of a phase transition ||21|-p^. Bauer and Campi were the first to 
apply the methods used in percolation studies to analyze MF data [pT|-p5|. In percolation 
theory the moments of the cluster distribution contain the signature of critical behavior 

In recent years further progress was made by experiments in which practically all the 
fragments emitted in a given event were detected, thereby permitting complete reconstruction 



of MF events [HHSl. The ALADIN Collaboration studied the MF of 400A-1000A MeV 



Xe, Au, and U nuclei on various targets. Their results showed that fragment yields were 
independent of the entrance channel when the data were scaled for projectile or target mass 
lO). The EOS Collaboration studied the MF of lA GeV gold on carbon and analyzed the 



data using methods developed in the study of critical phenomena [|36H4^. Several critical 
exponents were determined and their values suggested that the MF of Au can be understood 
as due to a continuous phase transition ||43| , |4^ . The first results of the MF of lA GeV La 



and Kr on carbon have also been reported [^^ . The MF transition appears to involve the 



breakup of a nucleus to form several IMF's. The production of nucleons and light particles 
- the nuclear analog of a gas- occurs largely in the fragment deexcitation step [0. Along 
with the earlier inclusive studies, these experiments suggested that prior to MF the remnant 
formed in the prompt stage achieves thermal equilibrium. A two step process was proposed 
for the collision. This two step process is an idealization of a time dependent process. In 
the first prompt stage nucleons are knocked out of the participants. The emission of these 
prompt particle leaves an equilibrated remnant nucleus which undergoes deexcitation in a 
second step. 

In the EOS experiment prompt stage was separated from the MF stage by making a cut 
on the kinetic energy of light charged particles |3^. This separation was made possible by 
the ability of the EOS detector to measure nearly all the charged particles and fragments 
emitted in each event using reverse kinematics ||3^. Complete reconstruction and stage 
separation was therefore possible for a majority of the events. These reconstructed events 
were characterized by mass, charge, and excitation energy of the remnant. The resultant 
remnant then undergoes MF to produce excited fragments which then undergo secondary 
decay to produce the observed fragments The ALADIN experiment is also capable of 
measuring all the charged particles except for Z=l particles The ISiS Collaboration 

has also produced one of the most complete MF data sets in the high energy collision of p, p, 
and 7i~ with gold ||32|-^ . In the ISiS experiment the charge and mass of the excited source 



were obtained on an event by event basis by subtracting the non equilibrium particles from 
the target charge and mass. However, the non-equilibrium particles could only be separated 



from the thermal particles by means of a parametrization involving two component moving 
source fits. 



The EOS ||3^ data can be compared with statistical models. There are several statistical 



models which have been used to study MF, [ll|,|3|,|^j48H56| but the most widely used are 
the Statistical Multifragmentation Model (SMM) [p|,^8|-pO|l and Monte Carlo microcanonical 



model (MMMC) [0,0. The theoretical interest in MF is not confined to only statistical 
models but several other approaches to study MF have been carried out e.g. percolation, 
lattice gauge and Ising models p2|j57| -|63| . A phenomenological droplet model ||6^, based on 
the Fisher droplet model, has also been used to describe the liquid-gas phase transition 



in nuclear reactions. Here we shall compare the experimental data with the predictions of 
the SMM model in the manner presented in previous publications [^^|65| , |66| . We found 



that this model is in good agreement with a variety of results for the MF of lA GeV Au 
p7| | . SMM requires mass, charge and excitation energy of the remnant as input. Using the 
experimental remnants from Au+C, the parameters of the SMM were fixed based on the 
agreement between SMM and data. Some comparisons with the data for La+C and Kr+C 
were also made in later work using the same set of parameters as in case of Au+C [^. 

This paper deals with the MF of lA GeV Au, La, and Kr on carbon. Section II gives 
a brief summary of the experimental remnant properties. The fragment properties are dis- 
cussed in Section III. Section IV describes SMM and gives a comparison with the data. 
The search for the critical transition is discussed in Section V. The determination of various 
critical exponents is given in Section VI. Energy fluctuations and heat capacity analyses are 
discussed in Section VII. Section VIII dwells on the nature of the phase transition in Au, La 
and Kr. Conclusions are given in Section IX. 

II. PROPERTIES OF THE REMNANT 

The reverse kinematic EOS experiment was performed with lA GeV ^^^Au, ^^^La, and 
^^Kr beams on carbon targets. The details of the experiment are given in our earlier publica- 



tions [^,0. Only a brief description will be given here. The experiment was done with the 
EOS Time Projection Chamber (TPC) p7| , p5| and a multiple sampling ionization chamber 
(MUSIC II) ||69|. The TPC provided almost 4 vr solid angle coverage in the center-of-mass 
system. Three-dimensional tracking and charged particle identification permitted momen- 
tum and energy reconstruction of fragments with charges in the range of 1 < Z < 8. Particle 
identification was based on specific energy loss along particle tracks. MUSIC II detected and 
tracked fragments with charges 8 < Z < Z^gam- The excellent charge resolution of this 
detector permitted identification of all detected fragments. 

The analysis presented here is based on fully reconstructed MF events for which the total 
charge of the system was taken as 79 < Z < 83, 54 < Z < 60, and 33 < Z < 39 for Au, 
La, and Kr, respectively p9|j45|] . Approximately ~ 32000, 26000, and 42000 events met the 
above criteria for Au, La and Kr, respectively. 
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FIG. 1. Average remnant mass, charge, and total charged-particle multiplicity as a function of 
thermal excitation energy E^f^ from the MF of Au, La and Kr. 

The remnant refers to the equilibrated nucleus formed after the emission of prompt 
particles. This remnant then undergoes MF. The charge and mass of the remnant were 
obtained by removing for each event the total charge of the prompt particles. In order to 
obtain the mass, the number of prompt neutrons was estimated by means of the ISABEL 
cascade calculation |70,7l|. The excitation energy per nucleon of the remnant, E*, was based 
on an energy balance between the excited remnant and the final stage of the fragments |j7^ 
for each event, as discussed in detail elsewhere P^,HH. 




In our previous work we have shown that some of the excitation energy of the remnant in- 
cludes a nonthermal component which may be ascribable to expansion energy, [P^ ,iT, 



The thermal excitation energy i?^*^ of the remnants is obtained as the difference between E* 
and Ex- This energy is an important quantity both for input to SMM and for the physics 



analysis of the data. Fig.[^ show plots of Aremn, Zremn and total charged particle multiplicity, 
m, for Au, La, and Kr, respectively, as a function of E^f^ • 

The plots in Fig.|l] show the nature of the remnant originating from the different system 
sizes. In case of Au and La one finds that the remnant size at the highest multiplicity is 
nearly 50% of the initial remnant size, while for Kr the remnant loses up to ~ 70 % of the 
initial mass. The range of E^f^ is also different in the three cases. For Kr the remnant reaches 
E^f^ as high as ~ 20 MeV/nucleon, while for Au and La the maximum E^f^ are ~ 9 and ~ 
12 MeV/nucleon, respectively. Figure shows the variation of the average m with E^i^, but 
does not give an indication of event-to-event fluctuations. Fig.^ shows a contour plot of m 
vs -E^. These two quantities are fairly closely correlated for Au and La and confirm the 
linear variation shown in Fig.|l|. The correlation is much poorer for Kr and a broad range of 
energies corresponds to a narrow range of m values above 8 MeV/nucleon. 



III. FRAGMENT PROPERTIES 

We first examine the second stage fragment multiplicities, mg. Fig.|^ shows a plot of mg 
vs E^f^ for Au, La, and, Kr. The mg for Au and La increase linearly with E^f^, while for Kr mg 
increases slowly up to E^f^r^ 8 MeV/nucleon and remains essentially constant at higher E^j^. 
The constancy of above 8 MeV/nucleon suggests that the system may be disintegrating 
into individual nucleons and light particles, which suggests that the vaporization process 
becomes dominant. This interpretation is confirmed by the variation with E^f^ of the size of 
the largest fragment, Amax, in each event as shown in Fig.^. It is observed for Kr that above 
~ 8-9 MeV/nucleon the value of Amax is nearly constant with a value of < 6. 
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FIG. 2. m as a function of E^f^ for Au, La, and, Kr. There are ten equidistant contours. 

The distribution of IMFs is shown in Fig.^ for all three systems as a function of energy. 
In order to exclude fission and/or the largest fragment for all three projctiles we define IMFs 
as having nuclear charges ranging from Z=3 to Z=Zproj/4:. Fig.|] shows that for all the three 
systems the peak in IMFs is ~ 8 MeV/nucleon although the IMF peak for Au is not well 
defined as i?^*^ > 8 MeV/nucleon was hardly achieved for this projectile. The drop in IMF 
yields for Au has been seen when the data were plotted as a function of m rather than E^f^ 
because of the dispersion in m values at a given E^f^ |^5| . Note that the number of IMFs from 
Kr drops to substantially less than 1 above 10 MeV/nucleon. This is another confirmation 
of the results given in the preceding section showing the vaporization of the Kr remnants at 
high E*f^. 
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FIG. 3. 1712 as a function of E"^^ for Au, La, and, Kr. 

IMF production has been studied as a function of projectile bombarding energy by the 
ALADIN Collaboration [p^ -|30|. They showed that the excitation energy dependence of the 
average IMF number for Xe, Au, and U projectiles when scaled by charge of the emitting 
source is similar, suggesting that the IMF production mechanism is independent of the 
entrance channel. A universal scaling for IMF production has also been seen from a wide 
range of source masses (35 -190 nucleons) produced in reactions with energies from 35 to 
600 MeV/nucleon 
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FIG. 4. Ajnax as a function of £^^j.^ for Au, La, and, Kr. 



IV. THE SMM MODEL AND COMPARISON WITH DATA 

The variable volume version of SMM has been used in this work for comparison with the 
data. This variable volume corresponds approximately to the condition of constant ( close 



to zero) pressure [0. SMM is a statistical description of the simultaneous breakup of an 
expanded excited nucleus into nucleons and hot fragments Individual fragments at 

normal nuclear density are described with a charged liquid drop parameterization. The free 
energy of a fragment is used to determine the fragment formation probability. This solution 
explicitly assumes the inhomogeneous nature of the hot MF final state. 
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FIG. 5. Average n(IMF) as a function of E^f^ for Au, La, and, Kr. 

Light fragments with Z < 3 may also be present in the hot MF final state. For the 
Z > 3 fragments, a quantum mechanical description is used for the temperature dependent 
volume, surface, and translational free energy of the fragments. The temperature indepen- 
dent parameters are based on the coefficients of the semiempirical mass formula. The critical 
temperature, at which the surface tension of neutral nuclear matter droplets would go to 
zero, is in the range suggested by infinite neutral nuclear matter calculations [ [f5[] . 

The two important parameters of the model are the Coulomb reduction parameter, k, and 
the inverse level density parameter eo |^^. The k parameter was fixed with the comparison 
of the measured free volume from SMM to that of the initial volume obtained in the collision 
process for the Au+C data. The details are given in ref. |^^. The only remaining parameter, 
eo, was obtained from the detailed comparison of SMM results with the various experimental 
fragment properties, e.g. second stage multiplicity, size of the biggest fragment, and IMFs 



using Au+C data |^^. In this work the same value of k=2 and eo=16 MeV has been used 
in SMM to compare results with data for La and Kr. The best agreement with the data was 
found by using the standard values of the parameters of the model. 
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FIG. 6. Second stage charged-particle multiplicity as a function of E"*^ for Kr, La, and, Au from 
data and SMM. 



The SMM results were also compared with the data on La and Kr with the same 
parameters as in Au, using reduced multiplicity (m/Zproj). Here we present a few compar- 
isons between data and SMM using E^f^. Fig.^ shows a plot of mg as a function of E^f^ for 
Au, La and Kr. The agreement between data and SMM is very good. The flatness of mg 
beyond E^f^r^ 8 MeV/nucleon observed for Kr is also seen in SMM. The size of the largest 
fragment from both data and SMM is shown in Fig.|^ for all the three systems. Very good 
agreement between the two is obtained over the entire range of E^f^. 
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FIG. 7. Size of the largest fragment as a function of E^f^ for Kr, La, and, Au from data and 



SMM. 



The average value of the total number of IMFs as a function of E^f^ is shown in Fig.|^ 
for both data and SMM. Both show the same initial increase in IMF production and a 
peak at approximately the same energy. SMM follows the trend in the Kr data in the 
vaporization regime, where the IMF multiplicity decreases to zero at the highest energy. 
SMM overestimates the number of IMF at the peak in all three cases. This difference is due 
to the fact that in SMM there is overproduction of Li and Be fragments at higher excitation 
energies. 
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FIG. 8. Average multiplicity of IMFs as a function of E^f^^ for Kr, La, and, Au from data and 
SMM. 

The above comparison between data and SMM was done with SMM calculations including 
the deexcitation of secondary fragments. SMM can also be stopped prior to secondary decay 
to obtain information about the SMM primary fragments. We refer to these results as 
SMMfiot while those following secondary decay are designated SMMcoid- As will be seen in 
the following sections SMMhot can be used to evaluate the nature of the phase transition. 
A full description of SMMhot for the remnant system with A=160 and Z=64 has been given 
in ref . ^ . 
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FIG. 9. Average multiplicity from SMM^ot and SMM^oid as a function of E^^^ for Kr, La, and. 



Au. 



The 7712 distributions as a function of E^^^ for both SMMcoid and SMMhot are shown in 
Fig.p. There is an increase in multiphcity for SMMcoid in all three cases. For SMM^ot ^2 
remains nearly unchanged up to a certain energy and then suddenly increases. The energy 
at which 7712 starts increasing corresponds to the MF threshold. The values of at which 
MF first occurs are ~ 4, ~ 5, and ~ 8 MeV/nucleon for Au, La and Kr, respectively. The 
threshold energy is lowest for the heaviest system because of the large Coulomb energy in 
the heavier remnant, which facilitates the breakup of the nucleus at a lower energy. There 
is a very narrow window for Kr between the multifragmentation threshold and the energy 
at which vaporization starts. This reduces the probability of IMF formation, as is evident 
from Fig.^ 



V. CRITICAL POINT DETERMINATION 



In our earlier publications we discussed the lA GeV Au on carbon data in terms of the 
theory of critical phenomena and several critical exponents were determined pB|,pB|p]|,^P 



We used the percolation technique applied to small lattices to study critical phenomena 
76| , [77|j35| . The method of moments analysis was used by several groups ||29| , |57| , |78| -^ to 



search for evidence of the liquid-gas phase transition in MF. Recently, we used this method in 
the analysis of La and Kr data |^ using charged particle multiplicity as the order parameter. 
In this work E*f^ is used as the order parameter. E*f^ is a more fundamental parameter than 
multiplicity for comparing the three systems. Here, we first use the combination of moments 
to find the signature of criticality in Au, La and Kr data. For example, the reduced variance 
72 = M2Mo/Mf, where Mi and M2 are the first and second moments of the mass distribution 
in an event and Mq is the total multiplicity including neutrons, is a useful quantity. The 
reduced variance 72 has a peak value of 2 for a pure exponential distribution, ua ~ e"""^, 
regardless of the value of a, but 72 > 2 for a power law distribution, ~ when r > 2 
and the system is large enough. Critical behavior requires that the peak value of 72 be larger 



than 2 21.221 
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FIG. 10. 72 as a function of Ef, for all the three systems from data and SMM, 



cold- 



The event by event analysis for 72 is shown in Fig.|TO| for all the three systems. The 
position of the maximum 72 value defines the critical point, i.e. the critical energy i?*, 
where the fluctuations in fragment sizes are the largest. The peak in 72 is well defined for 
La and Au. The value for E* is in good agreement with our earlier values, where was 
used as the order parameter and was obtained from E^^^ vs m plot m. In case of Kr 



Fig.|10| shows that there is no well defined peak in 72 and the distribution is very broad. A 



well defined peak is obtained in mc, but the value of 72 is always less than 2 ||46|| . 

Fig.0 also shows the 72 calculation using SMMcoid- The fission contribution to 72 has 
been removed both from the data analysis and SMM. The E* values obtained for data and 



SMMcoid are close to each other. There is a difference in the height of the peak in 72 for Au 
between data and SMMcoid- It is also important to note that at the peak 72 > 2 both from 
SMMcoid and data for Au and La and 72 < 2 for Kr. This suggests that the exponent r < 2 
for Kr as there is no enhancement of 72 in the critical region . One expects an enhancement 
of the moments in the critical region with r > 2 as in most critical phenomena [0,^ . 



In case of Au the 72 value remains above two for most of the excitation energy range. 
The E^f^ width over which 72 > 2 is smaller for La and disappears for Kr. The decrease in 



72 with decrease in system size, as observed in Fig. 10, is also seen in 3D percolation studies 



and the differences have been attributed to finite size effects |]83|j8^ . The results for E* are 
given in Table I along with the values of nic PBI . 
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It is of interest to see how the critical point determination from SMM^ot compares with 
that of SMMcoid- Since SMMcoid can reproduce the various features of the EOS data, one 
can use SMMhot to understand the MF mechanism. Fig. |Tl] shows a plot of 72 as function 
of Elfj^ from SMMcoid and SMM^ot- It is evident from the plots for Au and La that the 
peak in 72 occurs at the same i?^ for both SMMcoid and SMMhot- The distinct difference 
is in the height of the peak. For SMMhot the height of the 72 peak is smaller as compared 
to SMMcoid, but still above 2. This difference is mainly due to the increase in SMMcoid 



multiplicity as compared to SMMhot- In the case of Kr the 72 value is always less than 2 for 
both SMMcoid and SMMhot and the height of the 72 peak does not decrease from SMMcoid 
to SMMhot as for Au and La. A detailed calculation of critical exponents from SMM^ot has 
been given in our earlier publication These results clearly demonstrate that SMM^ot 



and SMMcoid behave in a similar way. The most important conclusion is that the values of 
the exponents do not change in going from SMMhot to SMMcoid- Thus, SMMhot results 
are consistent with a critical phase transition for Au and La and not for Kr. 

VI. CRITICAL EXPONENTS ANALYSIS 
A. T Exponent 

In the previous section it was shown, based on the 72 analysis, that Au and La should 
show a power law with exponent r > 2 for a continuous phase transition to be present, r 
can be obtained using the moments of the fragment mass distributions. 

Scaling theory relates the values of the critical exponents r and a to the moments 



Mfc of the mass distribution through 

Mfccx (1) 

where p is the bond breaking probability and at the critical point p=Pc in percolation. The 
values of r and a in the above equation are characteristic for the specific class of phase 
transition. For a transition of the 3D percolation type r = 2.2 and a = 0.45, while for a 
liquid-gas phase transition a value of r = 7/3 and a = 2/3. The 3D Ising values for r and 
a are 2.2 and 0.64 respectively. The r value is essentially the same for different universality 
classes. The above equation can be solved to get the value of r if the second (M2) and 
third moments (M3) of the fragment mass distributions are known. A plot of ln{M^) versus 
ln{M2) should give a straight line with a slope given by 

^ Mn{M,) ^ r-4 

A/n(M2) r-3 ^ ' 



Fig.0 shows a scatter plot of ln{M^) vs ln{M2) for all the three systems from data above 
E* [^. Also shown in Fig.|T2] are the results from SMMcoid- A linear fit to ln{M^) vs 
ln{M2) gives the value of r. The r values are shown in Table I both from data and SMM. A 
linear fit to the Kr data gives a value of r = 1.88 ± 0.08. This value is below the minimum 
value T > 2 expected for a continuous phase transition. This result for Kr is consistent with 



Fig.|TO|, which shows that 72 < 2 and hence r < 2 pl| , p2| . Similar results are also obtained 
for SMM^id. 
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FIG. 12. ln{M^) vs ln{M2) for Au, La and Kr above the critical energy. 
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FIG. 13. Power law fit to the fragment mass yield distribution at the critical point. The solid 



line corresponds to n^/ ~ qqA with qo=0.2 and r=2.2[43]. 
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TABLE L (Jrilicril pMrauit-lcrs [nmi cImIm Mud .S'il/il/,.„/,/ 
Parameter Audata Laoata Kraata Augmm Lusmm Krsmm Per" LG'' 



rric 


28±3 


24±3 


18±2 


26±3 


23±3 


17±2 






El 


4.5±0.5 


5.5±0.6 


6.5±1.0 


4.3±0.5 


5.3±0.6 


6.2ibL0 






T 


2.16±0.08 


2.10±0.06 


L88±0.08 


2.11±0.05 


2.05±0.05 


L81±0.06 


2.20 


2.21 


P 


0.32±0.02 


0.34±0.02 


0.53±0.05 


0.35±0.03 


0.37±0.03 


0.57±0.06 


0.44 


0.328 




0.22±0.03 


0.25±0.01 


0.50±0.01 


0.28±0.03 


0.29±0.05 


0.52±0.01 






7 


L4±0.3'^ 






1.02±0.23'^ 






1.76 


1.24 



1.32±0.15'' 1.20±0.08' 



a. Percolation, b. Liquid-Gas, c. reL[47], d. from P and P/^ ratio 



A one-parameter power law search also provides an alternative method for the determi- 
nation of r [Q. At the critical point the cluster distributions were described by n^/ ~ QqA^^ 



with r ~ 2.2 and go ~ 0.2, as seen in many universality classes Fig.|13| shows a plot of 



the fragment mass distribution at the critical energy for Au, La and Kr. The line corresponds 



to UAf = QqA with r=2.2 and qo=0.2 It is clear from the plot that Kr data do not 



follow the above power law with the critical values of go aiid r. 

B. /3 exponent 

Using the methods developed in percolation studies the value of the exponent /? can be 
obtained for the MF of lA GeV Au+C |3^. It is related to the size of the largest cluster by 
the relation 

Ama. ~ \ef (3) 

where e = p — Pc and e > 0. In the MF case p and Pc have been replaced by and E*. In 
the infinite lattice, the infinite cluster exists only on the liquid side of Pc- In a finite lattice 
a largest cluster is present on both sides of the critical point, but the above equation holds 
only on the liquid side. When e < , no infinite cluster exists and the size of the finite 
cluster is given by 

A ~ |e|-(/3+7) (4) 
where 7 is another critical exponent and related to the second moment, M2 

M2 ~ |e|-^ (5) 



Fig.0 shows a plot of In(Amax) vs — E* from Au, La and Kr. The values of (3 
obtained from the fit are given in Table I. The values of (3 for Au and La are 0.32±0.02 
and 0.34±0.02, respectively, and close to the value of 0.33 predicted for a liquid-gas phase 
transition. The value of f3 for Au is in agreement with our earlier reported value [^. In 



case of La (3 is close to the value obtained for Au, while the value of (5= 0.53±0.05 for Kr is 
much higher than that of Au and La. 

It is interesting to note that the SMM calculation also gives the same results as obtained 



in data. Fig.0 shows results from SMM as open circles for all the three systems. These 

comparisons with SMM are important, as this will help us to probe the order of phase 

transition using SMM in the three experimentally studied systems. 
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FIG. 14. ln{Amax) vs ln\E^f^ — E*\ for Au, La and Kr below the critical energy for exponent /5 
determination. 



Campi pT| , ^ also suggested that the correlation between the size of biggest fragment. 



Amax, and the moments in each event can measure the critical behavior in nuclei. Fig. 15 



shows a plot of the logarithm of A^ax in each event versus the logarithm of the second 
moment M2 for all the three systems. This plot is generally called a Campi scatter plot and 
has been successfully used in many studies |]85|j86[] . The two branches corresponding to the 
under-critical ( upper branch ) and over-critical (lower branch) events are clearly seen for 



Au and La. The Amax - M2 correlation is quite broad for Kr and fills most of the available 
phase space. The lower and upper branches seems to overlap and are not well separated. 
Studies on percolation lattices show similar behavior |8^. However, it is not possible from 
such a plot to locate the critical region in a precise and unambiguous manner. But, if we 
know the critical point from some other method, then the Amax- M2 correlation can be used 
to calculate the ratio of critical exponents /?/7 from the slope of the upper branch. A similar 



behavior is obtained from SMM as shown in Fig. ^ for all the three systems. 
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Table 1 shows the jS/'-f values for Au and La from the linear fit to the upper branch of 
Fig.|T5|. For the fitting purpose an average value of ln{Amax) was obtain for each value of 
/n(M2). 

We have mentioned only a few exponents in this analysis and compared them with SMM. 
The aim is to show that Kr is different than Au and La. The two exponents r and (3 serve 
this purpose very well. Knowing (3 and jS/j, '-f can be obtained. Results are listed in Table 
I and the value for Au agrees with the published value. Table I also gives the values of 
the critical exponents for the percolation and liquid-gas phase transition in 3D systems. No 



attempt has been made in this paper to independently determine other exponents hke 7 and 
a. All the exponents for Au and their comparison with SMM have been reported in various 
pubhcations |6|,|8|,|0|,|7 . 



VII. ENERGY FLUCTUATIONS AND HEAT CAPACITY ANALYSIS 

In a recent study experimental evidence of a liquid-gas phase transition in MF was offered 
by analyzing the fluctuations in the total fragment energy |9^j9^. It has been shown that 
for a given total energy the average partial energy stored in a part of the microcanonical 
system is a good thermometer while the fluctuations associated with the partial energy can 
be used to determine the heat capacity, which is negative for a first order transition and 
positive for a second order transition. 
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FIG. 16. With no selection on the remnant mass: a). Partial energy Ei as a function of E* b). 
Fluctuation in Ei c). The variance (a/Ei)^ and d). Heat capacity per nucleon for Au . 



A negative heat capacity was obtained in the reaction Au+Au at 35A MeV, providing the 
direct evidence of a first order liquid-gas phase transition [Q. We have analyzed our data 
for Au and Kr using energy fluctuations. From the total excitation energy E* the Coulomb 



and expansion energy components were removed to obtain energy Ei |]39| , [45 . 
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FIG. 17. With a tight cut on remnant mass: a). Partial energy Ei as a function of E* b). 
Fluctuation in Ei c). The variance {a/EiY and d). Heat capacity per nucleon for Au. 

The fluctuations in Ei were studied as a function of E* . Figs. |T^(a) and (b) show the 
spread in Ei and its standard deviation a respectively, as a function of E* . The reduced 



variance (a/Ei)"^ is shown in Fig. |T^(c). The heat capacity Ct is given by 



a = C,V[Ci - (a/T)'] 



(6) 



where Ci is the canonical specific heat. It was argued that a negative Ct in the critical region 
is a signature of a first order phase transition [92,93]|. In Fig. 0(d) Ct is shown as a function 
of E*. There is a negative heat capacity at E* ~ 2-3 MeV/nucleon. 

However, a different result is obtained if there is a selection on the remnant mass. 
Figs. |TB|(a) and (b) have no remnant mass cut and the average mass is 170 ± 19. This 
causes large fluctuations in Ei. Figs. |T^(a)-(d) show similar plots as in Figs. |TB|(a)-(d) with 
the remnant restricted to 161 ± 3. Two important results come out of these plots. First, 
a is much smaller when the remnant mass cut is applied. Second, Ct is always positive. 
Our analysis for Au data shows no negative heat capacity. Based on our statistical and 
thermodynamic analysis, we see indications for a continuous phase transition in Au data 
ll36|,|38|M0|,|47ll . The results for La are similar to those shown above for Au. 



A similar analysis was performed on Kr data and results are shown in Figs. |l8|(a)-(d) 
and Figs, ^(a)-(d) for remnant masses 58 ± 14 and 56 ± 3 respectively. We do not see the 
signature of a first order transition even in Kr if the remnant mass is selected in a narrow 
bin. On the other hand a negative value of Ct is obtained if there is no selection on the 
remnant mass. 

Here, we emphasize again that the event-by-event complete reconstruction of remnant 
mass and excitation energy is very important for the above analysis. This is possible only in 



high energy reverse kinematic asymmetric collisions |3^. In case of the EOS data, though 
there is a positive specific heat, there are no large scale fluctuations in it to define the critical 
energy. Only Au has a peak in Ct-, which is around 2-3 MeV/A. There is no peak in Ct for 
La and Kr. Thus the analysis of EOS data for all three systems based on heat capacity is 
questionable. 
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FIG. 18. With no selection on the remnant mass: a). Partial energy Ei as a function of E* b). 
Fluctuation in Ei c). The variance {a/Ei)"^ and d). Heat capacity per nucleon for Kr 
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VIII. NATURE OF PHASE TRANSITION 



In the previous sections we have shown that the MF of Kr is different than that of Au 
and La. Based on the experimental resuhs it is not possible to decide the order of phase 
transition in Kr. The statistical analysis suggests a continuous phase transition for Au and 
La. 
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In our earlier publication |^ and in this paper as well we have shown that SMM can 
reproduce the various features of EOS data, including the critical exponents. A single set of 
parameters of the model was used to describe the data for all the three systems, Au, La, and 
Kr. The nature of the phase transition in SMM was analyzed using SMMhot, i-e. fragments 
formed before secondary decay [^. The microcanonical temperature was obtained for the 
MF system and a caloric curve was constructed p|, ^p4| for A=160, 130, 100, and A=70. A 
backbending in the caloric curve is an indication of a first order phase transition and leads 
to a negative specific heat. For A=160 there is a positive peak in the specific heat vs energy 
plot, consistent with our energy fluctuation analysis on our data, using a mass cut. 
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FIG. 20. P = l/Tfyreak a function of E* from SMM^ot using experimental remnant as input 
to SMM. The graphs are shown for different mass selection on the remnants from MF of Au. 



The back-bending in the caloric curve is sensitive to the remnant mass. We have used 
the remnant distributions from Au+C and Kr+C data to study the remnant mass effect 
on the caloric curve. Figs.pO|(a)-(d) show a plot of the reciprocal of the SMM^ot breakup 
temperature (3 = 1/Threak, for Au. A distinct pattern is evident. A small back bending starts 
appearing at the very low mass cut. For heavier masses there is no back-bending. Since the 
remnant mass distribution from Au+C interaction is dominated by heavier remnants, the 
effect of lighter remnants on back-bending is almost negligible. Fig.|2^ suggests that, accord- 
ing to SMM, the MF transition in Au changes from first order to second order depending 
on the remnant mass. For Kr, the back-bending is very sensitive to the particular mass 
cut. If there is no selection on mass cut the caloric curve is different than those with mass 



cuts as shown in Figs. |2T|(a)-(d). However, back-bending is always present. Thus, based on 
the above reasoning we can say that in case of Kr there is no continuous phase transition. 
Rather, analysis of the caloric curve argues for a first order phase transition in Kr. 
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FIG. 21. (3 = l/Tfej-eafc as a function of E* from SMMj^^t using experimental remnant as input 
to SMM. The graphs are shown for different mass selection on the remnants from MF of Kr. 



So far we have shown that the Kr data are different than the Au and La data. The 
temperature of the MF system provides further insight. One can obtain the freeze out 
temperature of the cold fragments in data using the ratio of light fragment isotopic yields 
Il39| , |45|j96| . In Fig.^(a) the isotope ratio temperature (Th^e-dt), as obtained from ^H/^H to 
^He/'^He yield ratio at the critical point E*, is shown for Au, La and Kr systems as function 
of the linear size of the system (L= A^J^^). 
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FIG. 22. a). Tue-DT temperature at -E* as a function of system size from data b). Critical 

temperature from Ising model calculation and c). Critical bond formation probability plotted as 
1/pc for different lattice sizes. 

Thc-dt decreases with an increase in system size. Recently, in an another study it has 
also been shown that critical temperatures and excitation energies decrease with increasing 
This result can be attributed to the higher Coulomb energy in Au as 



system size 



compared to Kr, which in turn shifts the MF transition to lower temperature. The behavior 
of Thb-dt with system size is different than those observed either in percolation or Ising 
model studies. Fig. ^(b) shows a plot of the critical temperature from a 3D Ising-type model 
with fixed density for different lattice sizes. In the Ising calculation the trend is different 
than the one observed in Fig. p^(a). Thus, for neutral matter the critical temperature 
increases with an increase in system size. We can also compare the results from data and the 
Ising model calculation with percolation studies. In percolation the critical probability in 



bond building percolation decreases with increase in the system size [[77|,p7|. Fig. |22|(c) shows 



1/ Pc percolation lattice size, where Pc is inversely proportional to Tc ||9^. Thus, MF of 
Au, La and Kr is different than 3D Ising and percolation models. As mentioned earlier, the 
finite size affects only non-thermodynamical quantities e.g. 72, A^ax, etc. and not energy or 
temperature. 

The critical energy is shown in Fig.|23|( function of system size. The decrease in 

E* with increase in system size is evident from the figure. SMM also predicts the critical 
energy of the MF transition and this energy is shown in Fig.^(a) along with the data. This 



energy is obtained at the hot fragment stage from the peak in 72 as shown in Fig.|Tl]. Both 
data and SMM E* are in good agreement. The breakup temperature Tsmm, as obtained 
from SMMhot, along with the Thc-dt from data is shown in Fig.^(b). There is a decrease 
in both temperatures with increase in system size. It is apparent that Thc-dt is about 1 
MeV lower than the SMM temperature. This difference is due to the fact that Thb-dt is 
measured after secondary decay has taken place, while Tsmm corresponds to the breakup 
configuration. Note that Thc-dt tracks Tsmm with system size at the critical point. SMM 
indicates that the decrease in both Tsmm and E* with increasing system size is due to the 
increase of Coulomb energy. This result suggests that Coulomb energy plays an important 
role in the MF of nuclei. It is interesting to note that the critical energy obtained both from 
SMM and data corresponds to the MF threshold as shown in Fig.^. The significance of this 
result remains to be determined. 

The effects of finite size and Coulomb force on the MF have been studied by several 
workers and it is found that the decrease in critical temperature with increase in system size 



is primarily due to the Coulomb energy 1 98 -IOC]. The result from one such calculation is 
shown Fig.^(b) as Tumu ■ 
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The microcanonical Metropohs Monte Carlo (MMMC) calculations have emphasized 
that MF is controlled by the competition between long range Coulomb forces and finite size 
effects ( especially surface energy). Finite size effects in models with only short range forces 
predict an increase in the critical temperature as the system size increases, as is evident from 
percolation |^ and Ising model studies |^ (see Fig.^2|). Since the experimental temperature 
exhibits the opposite dependence on system size, it is apparent that Coulomb effects are more 
important than finite size effects. For finite neutral matter the critical temperature (Tc) is 
expected to be ~ 15-20 MeV HllOl. The observed for A=160 is ~ 6 MeV. Compared 
to finite uncharged nuclei, the presence of Coulomb energy plays a role in lowering the 
excitation energy needed to reach the regime where critical signatures are observed. In the 
smaller Kr system there is less Coulomb energy in the initial remnant state. Achieving 



multifragmentation in this system requires greater excitation energy /nucleon compared to 
Au and La (as shown in Fig. |2^(a)) and as a result, the dynamics of the ensuing disassembly 
may not take this system near its critical regime. 
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FIG. 24. Ml distribution for SMMhot for Au, La and Kr. 

It has also been suggested on the basis of a caloric curve that MF of Au is a first order 
phase transition [|102|| . This observation was based on the the fact that the temperature 
remains constant as energy is increased between 3 and 10 MeV/nucleon. The temperature 
again starts increasing beyond 10 MeV/nucleon as a gas phase consisting of a mixture of 
nucleons and few light particles is created. According to SMM there are few nucleons in 
the hot system. Fig.^ shows a plot of first moments of the fragment yield distribution as a 
function of E^,^ for Au, La and Kr. The figure shows that even for E^f^ > 10 MeV/nucleon 
very few particles with A < A are produced. Most of the remnant mass is in fragments with 
A > 4. This argues against the coexistence of the constant density liquid and gas phases 



in SMM. The cooling of the hot fragments produces a large number of final state nucleons 
and light particles as shown in Fig.EH The IMFs survive the cooling process at least for Au 
and La and identify the MF transition. In case of Kr the excitation energies of IMF from 
SMMhot are very high and very few are seen in the cold stage. Thus the MF signal could 
be washed out. 
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FIG. 25. Ml distribution for SMMcoid for Au, La and Kr. 



IX. SUMMARY AND CONCLUSIONS 

In the present work we presented and analyzed the data from the MF of lA GeV Au, La 
and Kr on carbon. The mass, charge and excitation energy of the remnant were determined 
in each event. The thermal excitation energy was obtained after the expansion energy was 
determined on the basis of energy balance. The multiplicity distribution from MF of Kr 
shows a saturation beyond E^f^ of 8 MeV/nucleon, indicating that the vaporization process 



has started. 

A comparison of the data with the variable volume version of SMM, using the same 
values of parameters of the model for Au, La and Kr show a very good agreement with 
various distributions e. g. charged particle multiplicity, size of biggest fragment, number of 
IMFs, etc. The power law behavior seen in data for Au and La with r > 2 is also observed 
in SMM when cold fragments are analyzed. The hot fragment analysis from SMMhd for 72 
gives a value of E* which is in agreement with the value of E* obtained from SMMcoid- 

The Au and La data give indications of a continuous phase transition. A first order 
phase transition has been predicted for Kr using SMM. However, the data cannot be used 
to distinguish between the two because of finite size effects. 

The SMM calculation offers the best way to determine the the nature of phase transition. 
A back bending in temperature vs energy plot is a sign of negative specific heat and hence 
a first order phase transition. Such a result is found for Kr but not for Au, except for the 
lightest remnants, the MF contribution of which is negligible. 

The temperature obtained from the isotope ratio analysis decreases with an increase in 
system size at the critical point. The break-up temperature obtained from SMM also follows 
the same trend, but higher by ~ 1 McV/nuclcon. The decrease in temperature with increase 
in system size is an important result as it is in the opposite direction from what is observed in 
either in 3D percolation or Ising model studies for finite size neutral matter. This shows that 
if finite size were the only effect in MF then the correlation between temperature and system 
size would have been the same in data as in percolation and Ising model studies. Thus, the 
long range Coulomb force is shown to be the dominant factor in MF. This conclusion is also 
supported by Hartree-Fock calculations. 

In conclusion, this is the first work in which the nature of the phase transition in MF has 
been explored using three systems of different size. The experimental results in conjunction 
with SMM provide the order of the phase transition in Au, La and Kr. The values of critical 
exponents r, /3 and 7, which are close to the values for liquid-gas system, along with nearly 
zero latent heat suggest a continuous phase transition in Au and La. The back bending in 



the caloric curve for Kr suggests the presence of latent heat, which is consistent with a first 
order phase transition. We emphasize again here the important role played by the Coulomb 
energy. The Coulomb expansion energy reduces or eliminates the latent heat and changes 
the nature of the phase transition. 

This work was supported by the U. S. Department of Energy. 



orv 



REFERENCES 

[1] D. H. E. Gross ,Rep. Prog. Phys. 53, 605 (1990). 

[2] L. G. Moretto and G. J. Wozniak 43, 379 (1993). 

[3] J. Bondorf et al, Phys. Rep. 257, 133 (1995). 

[4] D. H. E. Gross , Phys. Rep. 279, 119 (1997). 

[5] J. Richert and P. Wagner, Phys. Rep. 350, 1 (2001). 

[6] A. Bonasera et al, La Rivista del Nuovo Cimento, della Societ a Itahana di Fisica(2000). 

[7] S. Das Gupta, A. Z. Mekjian, and M. B. Tsang , |nucl-th/ 0009033 . 

[8] J. E. Finn et al, Phys. Rev. Lett. 49, 1321 (1982). 

[9] R. W. Minich et al, Phys. Lett. B118, 458 (1982). 
[10] M. E. Fisher, Physics 3, 255 (1967). 
[11] P. J. Siemens, Nature 305, 410 (1983). 

[12] T. J. Schalgel and V. R. Pandharipande, Phys. Rev. G 36, 162 (1987). 

[13] H. MuUer, and B. D. Serot, Phys. Rev. G 52, 2072 (1995). 

[14] K. G. Ghase and Z. Mekjian, Phys. Rev. Lett. 75, 4732 (1995). 

[15] A. M. Poskanzer et al, Phys. Rev. G 3, 882 (1971). 

[16] G. D. Westfall et al, Phys. Rev. G 17, 1368 (1978). 

[17] A. S. Hirsch et al, Phys. Rev. G 29, 508 (1984). 

[18] M. Mahi et al, Phys. Rev. Lett. 60, 1936 (1988). 

[19] N. T. Porile et al, Phys. Rev. G 39, 1914 (1989). 

[20] G. J. Waddintgton and P. S. Freier , Phys. Rev. G 31, 888 (1985). 



[21] X. Campi, J. Phys. A 19, L917 (1986). 
[22] X. Campi, Phys. Lett. B 208, 351 (1988). 

[23] X. Campi, Proceedings of the International School of Physics 'Enrico Fermi' Nuclear 
Collisions from the Mean-Field into the Fragmentation regime, CXII 331 (1991). 

[24] W. Bauer et al, Phys. Lett. B150, 53 (1985). 

[25] W. Bauer et al., Nucl. Phys. A452, 699 (1986). 

[26] D. Stauffer, Phys. Rep. 54, 1 (1979). 

[27] D. Stauffer and A. Aharony, Introduction to Percolation Theory ( Taylor and Francis 
, London, 1992). 

[28] J. Hubele et al, Z. Phys. A 340, 263 (1991). 

[29] P. Kreutz et al, Nucl. Phys. A 556, 672 (1993). 

[30] A. Schuttauf et al, Nucl. Phys. A 607, 457 (1996). 

[31] U. Jahnke et al., Phys. Rev. Lett. 83, 4959 (1999). 

[32] W.-C. Hsi et al, Phys. Rev. Lett. 79, 817 (1997). 

[33] L. Beaulieu et al, Phys. Lett. 163B, 159 (1999). 

[34] T. Lefort et al, Phys. Rev. C64, 064603 (2001). 

[35] M. Kleine Berkenbusch et al, Phys. Rev. Lett. 88, 022701 (2002). 

[36] M. L. Gilkes et al, Phys. Rev. Lett. 73, 1590 (1994). 

[37] J. A. Hauger et al, Phys. Rev. Lett. 77, 235 (1996). 

[38] J. B. Elliott et al, Phys. Lett. B381, 35 (1996). 

[39] J. A. Hauger et al, Phys. Rev. C57, 764 (1998). 



[40] J. B. Elliott et al, Phys. Lett. B418, 34 (1998). 

[41] J. Lauret et al, Phys. Rev. C57, R1051 (1998). 

[42] B. K. Srivastava et al., Phys. Rev. C60, 064606 (1999). 

[43] J. B. Elliott et al, Phys. Rev. C62, 064603 (2000). 

[44] J. B. Elliott et al, Phys. Rev. Lett. 85, 1194 (2000). 

[45] J. A. Hauger et al, Phys. Rev.C62, 024616 (2000). 

[46] B. K. Srivastava et al. Phys. Rev. C64, 041901 (R) (2001). 

[47] R. P. Scharenberg et al, Phys. Rev.C 64, 054602 (2001). 

[48] J. Bondorf et al, Nucl. Phys. A 443, 321 (1985). 

[49] J. Bondorf et al, Nucl. Phys. A 444, 460 (1985). 

[50] J. Bondorf et al, Nucl. Phys. A 475, 663 (1987). 

[51] W. A. Friedman, Phys. Rev. C42, 467 (1990). 

[52] D. Durand, Nucl. Phys. A 541, 266 (1992). 

[53] K. C. Chase and A. Z. Mekjian, Phys. Rev. C 49, 2164 (1994). 

[54] K. C. Chase, A. Z. Mekjian, and P. Bhattacharya, Phys. Rev. C 55, 1410 (1997). 

[55] S. J. Lee and A. Z. Mekjian, Phys. Rev. C 56, 2621 (1997). 

[56] S. Das Gupta and A. Z. Mekjian, Phys. Rev. C 57, 1361 (1998). 

[57] W. Bauer, Phys. Rev. CSS, 1297 (1988). 

[58] J. Pan and S. Das Gupta, Phys. Rev. C 51, 1384 (1995). 

[59] X. Campi and H. Krivine, Nucl. Phys. A 620, 46 (1997). 

[60] S. Das Gupta et al, Nucl. Phys. A 621, 897 (1997). 



[61] J. Pan and S. Dasgupta, Phys. Rev. Lett. 80, 1182 (1998). 
[62] J. M. Carmona et al, Nucl. Phys. A 643, 115 (1998). 
[63] F. Gulminelli and Ph. Chomaz, Phys. Rev. Lett. 82, 1402 (1999). 
[64] A. L. Goodman et al, Phys. Rev. C 30, 851 (1984). 

[65] R. P. Scharenberg, Proceedings of the International Workshop XXVII on Gross Prop- 
erties of Nuclei and Nuclear Excitations, Hirschegg, Austria, p 237 , GSI 1999. 

[66] B. K. Srivastava, Proceedings of the International Workshop XXVII on Gross Proper- 
ties of Nuclei and Nuclear Excitations, Hirschegg, Austria, p 247 , GSI 1999. 

[67] G. Rai et al., IEEE Trans. Nucl. Sci. 37, 56 (1990). 

[68] H. Wieman et al, Nucl. Phys. A 525, 617c (1991). 

[69] G. Bauer et al, Nucl.Instr. Metds. Phys. Res. A 386, 249 (1997). 

[70] Y. Yariv and Z. Fraenkel, Phys. Rev. C20, 2227 (1979). 

[71] Y. Yariv and Z. Fraenkel, Phys. Rev. C24, 488 (1981). 

[72] D. Cussol et al, Nucl. Phys. A 561, 298 (1993). 

[73] L. Beaulieu et al, Phys. Rev.C54, R973 (1996). 

[74] I. N. Mishustin, Nucl. Phys. A 630, 111c (1998). 

[75] D. G. Ravenhall et al, Nucl. Phys. A 407, 572 (1983). 

[76] J. B. Elliott et al, Phys. Rev. C49, 3185 (1994). 

[77] J. B. Elliott et al, Phys. Rev. C55, 1319 (1997). 

[78] Bao-An Li et al, Phys. Lett. B303, 225 (1993). 

[79] H. W. Barz et al, Nucl. Phys. A 561, 466 (1993). 



[80] M. Belkacem et al, Phys. Rev. C54, 2435 (1996). 

[81] P. F. Mastinu et al, Phys. Rev. C57, 831, (1998). 

[82] M. D'Agostino et al, Nucl. Phys. A 650, 329 (1999). 

[83] X. Campi and H. Krivine, Nucl. Phys. A 545, 161c (1992). 

[84] X. Campi and H. Krivine, Z. Phys. A344, 81 (1992). 

[85] H. R. Jaqaman and D. H. E. Gross, Nucl. Phys. A 524, 321 (1991). 

[86] M. Belkacem et al Phys. Rev. C52, 271 (1995). 

[87] H. E. Stanley, Introduction to phase transition and critical phenomena, (Clarendon, 
Oxford, 1971). 

[88] M. N. Barber, Phase tensitions and critical phenomena, Vol.8 , Eds. C. Domb and J. 
L. Lebowitz ( Academic press Inc., London, 1983). 

[89] G. Grimmett, Percolation, (Springer, New York, 1989). 

[90] J. M. Yeomans, Statistical mechanics of phase transitions, (Clarendo, Oxford, 1993). 

[91] X. Campi, Nucl. Phys. A495, 259c (1989). 

[92] Ph. Chomaz and F. GulmineUi, Nucl. Phys. A 647, 153 (1999). 

[93] M. D'Agostino et al, Phys. Lett. B 473, 219 (2000). 

[94] A. Huller , Z. Phys. B 93, 401 (1994). 

[95] J. B. Natowitz et al , |nucl-ex/ 0106016. 



[96] S. Albergo et al, Nuovo Cimento A89, 1 (1985). 

[97] D. W. Heerman and D. Stauffer , Z. Phys.B44, 339 (1981) 

[98] H. R. Jaqaman et al, Phys. Rev. C29, 2067 (1984). 



[99] P. Bonche et al, Nucl. Phys. A436, 265 (1985). 
[100] S. Levit and P. Bonche , Nucl. Phys. A437, 426 (1985). 
[101] J. M. Lattimer et al, Nucl. Phys. A432, 646 (1985). 
[102] J. Pochodzalla et al., Phys. Rev. Lett. 75, 1040 (1995). 



A cr 



